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ABSTRACT 

Numerical simulations of the interiors of radiation dominated accretion disks show that significant density 
inhomogeneities can be generated in the gas. Here, we present the first results of our study on X-ray reflection 
spectra from such heterogeneous density structures. We consider two cases: first, we produce a number of toy 
models where a sharp increase or decrease in density of variable width is placed at different depths in a uniform 
slab. Comparing the resulting reflection spectra to those from an unaltered slab shows that the inhomogeneity 
can affect the emission features, in particular the Fe Ka and O VIII Lya lines. The magnitude of any differences 
depends on both the parameters of the density change and the ionizing power of the illuminating radiation, but 
the inhomogeneity is required to be within ~ 2 Thomson depths of the surface to cause an effect. However, only 
relatively small variations in density (on the order of a few) are necessary for significant changes in the reflection 
features to be possible. Our second test was to compute reflection spectra from the density structure predicted 
by a simulation of the non-linear outcome of the photon bubble instability. The resulting spectra also exhibited 
differences from the constant density models, caused primarily by a strong 6.7 keV iron line. Nevertheless, 
constant density models can provide a good fit to simulated spectra, albeit with a low reflection fraction, between 
2 and 10 keV. Below 2 keV, differences in the predicted soft X-ray line emission result in very poor fits with a 
constant density ionized disk model. The results indicate that density inhomogeneities may further complicate the 
relationship between the Fe Ka equivalent width and the X-ray continuum. Calculations are still needed to verify 
that density variations of sufficient magnitude will occur within a few Thomson depths of the disk photosphere. 

Subject headings: accretion, accretion disks — instabilities — line: formation — radiative transfer — X-rays: 
general 



1. INTRODUCTION 

The reprocessing of hard X-rays in relatively cold, optically 
thick material is a major component of the current X-ray phe- 
nomenology of active galactic nuclei (AGN) and Galactic black 
hole candidates (GBHCs). The evidence for this is based on the 
shape of the observed spectral continuum. Starting with Ginga 
observations in the late 1980s and continuing on with data col- 
lected from ASCA, RXTE and BeppoSAX, it was found that the 
continuum of almost all type 1 Seyfert galaxies exhibit a spec- 
tral hardening beyond ^10 keV (Pounds et al. 1990; Nandra & 
Pounds 1994) before rolling over at energies > 100 keV (Per- 
ola et al. 2002; Risaliti 2002). This fact combined with the 
near ubiquitous presence of a Fe Ka line (Nandra & Pounds 
1994; Reynolds 1997) is completely consistent with the predic- 
tions of Compton reflection of X-rays in dense, optically thick 
material (George & Fabian 1991; Matt, Perola & Piro 1991). 
Moreover, the equivalent width (EW) of the Fe Ka lines are 
usually large enough (EW ~ 50-300 eV) that the reflector must 
cover about half the sky as seen from the X-ray source (George 
& Fabian 1991), but remain out of the line of sight (in order 
to not produce significant absorption). In the complex environ- 
ment of the central engine of an AGN, there are many potential 
locations for the reprocessing gas, such as broad-line region 
clouds, and the obscuring material of the unification schemes. 
Indeed, reflection from one or both of these sites may be a com- 
mon occurrence in Seyfert galaxies (e.g. Yaqoob et al. 2001; 
Kaspi et al. 2001; Padmanabhan & Yaqoob 2003). However, 
in at least some objects (most famously, MCG-6-30-15; e.g., 
Fabian et al. 2002b) it was found that the Fe Ka line exhibits a 



broad red-wing, consistent with originating from the inner re- 
gions of the accretion disk (Fabian et al. 1989, 2000; Reynolds 
& Nowak 2003). Thus, in these sources, a detailed compari- 
son of the X-ray spectrum with models of Compton reflection 
may yield information on the structure, metallicity, and ion- 
ization state of the accretion disk at small Schwarzschild radii 
(Rs = 2GM/c 2 , where M is the black hole mass). Evidence 
for broad Fe Ka lines in GBHCs is becoming more common 
(e.g., Martocchia et al. 2002; Miller et al. 2002a,b,c, 2003), so 
the possibility of performing similar analyses also exist for the 
Galactic accreting black hole systems. 

Calculations of X-ray reflection from optically-thick material 
have increased in sophistication over the last decade. The initial 
Monte-Carlo calculations by George & Fabian (1991) and Matt 
et al. (1991) assumed a neutral, constant density slab, and were 
able to quantify the observed strength of the Fe Ka fluorescence 
line as a function of viewing angle and abundance. Soon after, 
Ross & Fabian (1993) and Zycki et al. (1994) allowed the gas 
to be ionized by the incoming X-rays and produced reflection 
spectra that included recombination lines and edges as well as 
fluorescence lines. Matt et al. (1993, 1996) studied the effects 
of ionization on the Fe Ka line, while Ross et al. (1999) clearly 
showed that the entire reflection continuum changed shape as 
a function of the ionization parameter (£ = 4irFx/nn, where 
Fx is the incident X-ray flux, and nn is the density of the re- 
processor). More recently, the assumption of constant density 
slabs has been replaced with more complicated density distri- 
butions, such as hydrostatic balance (Nayakshin, Kazanas & 
Kallman 2000; Nayakshin & Kallman 2001; Ballantyne, Ross 
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& Fabian 2001b; Rozanska et al. 2002) or constant pressure 
(Dumont et al. 2002). The major difference between these vari- 
able density models and the constant density ones is that under 
certain conditions (when the Compton temperature of the radi- 
ation field is very high) the illuminated atmosphere can be sub- 
ject to a thermal instability which can split the gas into an outer, 
low density, completely ionized zone and a deeper, high den- 
sity, completely recombined zone (cf., Krolik, McKee & Tarter 
1981). When this instability occurs the reflection spectrum ex- 
hibits only neutral features formed in the lower layer that are 
slightly broadened due to Compton scattering by the hot sur- 
face layer. Otherwise, the reflection spectra looks much like a 
diluted version of an ionized constant density model (Ballan- 
tyne et al. 2001b; Done & Nayakshin 2001). The features are 
weaker due to both Compton broadening by the hot surface gas, 
and the lower emissivity (proportional to density) of the emit- 
ting material. 

Over the last few years it has been possible to fit models 
which self-consistently predict both the reflection continuum 
and line emission to real X-ray data. The constant density mod- 
els of Ross & Fabian (1993) have been applied to both narrow- 
line and broad-line Seyfert 1 galaxies (Ballantyne, Iwasawa & 
Fabian 2001a; Orr et al. 2001; De Rosa et al. 2002; Longinotti 
et al. 2003) as well as GBHCs (Miller et al. 2003). While these 
are not as sophisticated as the hydrostatic models, they are rel- 
atively quick to compute and provide a good parameterization 
of the data which allows easy comparison between objects. On 
the whole, when combined with a relativistic blurring function, 
the models do provide a reasonable description of the hard X- 
ray spectrum for many accreting black holes. The exceptions 
to this, such as the narrow-line Seyfert 1 1H 0707-495 (Boiler 
et al. 2002) and MCG-6-30-15 (Ballantyne, Vaughan & Fabian 
2003) may indicate that more complicated geometries are re- 
quired than a simple flat accretion disk. 

As interesting as these fits are, a lingering question remains: 
are the constant density or hydrostatic assumptions represen- 
tative of real accretion disks? Where radiation pressure ex- 
ceeds gas pressure, as likely occurs in the inner regions of 
AGN disks accreting at more than 0.3% of the Eddington rate 
(Shakura & Sunyaev 1973), models with stress assumed pro- 
portional to total pressure are viscously (Lightman & Eard- 
ley 1974), thermally (Shakura & Sunyaev 1976), and convec- 
tively (Bisnovatyi-Kogan & Blinnikov 1977) unstable. These 
instabilities may lead to space or time variation in disk struc- 
ture. Furthermore, the stresses driving accretion in ionized, 
differentially-rotating disks are now thought to be due to mag- 
netic forces in turbulence resulting from the magneto-rotational 
instability or MRI (Balbus & Hawley 1991). In MHD simula- 
tions, gas motions in the turbulence vary over an orbital pe- 
riod (Balbus & Hawley 1998), which is similar to the timescale 
needed for establishing vertical hydrostatic balance (Frank, 
King & Raine 2002), Therefore, hydrostatic reflection models 
may not accurately represent the structure of turbulent accre- 
tion disks. Simulations of the turbulence which include ef- 
fects of radiation diffusion indicate that density fluctuations 
may exceed an order of magnitude (N.J. Turner et al. 2003). 
Radiation-supported disks with a magnetic field are subject to 
another dynamical instability, the photon bubble instability or 
PBI (Arons 1992; Gammie 1998). Photon bubble modes hav- 
ing wavelengths shorter than the gas pressure scale height can 
grow faster than the orbital frequency (Blaes & Socrates 2001, 
2003). The non-linear development of the instability may lead 



to trains of shocks propagating through the disk surface lay- 
ers (Begelman 2001), leading to a way of producing super- 
Eddington luminosities (Begelman 2002). Density contrasts 
between shocked and inter-shock regions can exceed one hun- 
dred (N.J. Turner et al. 2004, in preparation). Will density 
fluctuations caused by the MRI and PBI affect the reflection 
spectrum? Fabian et al. (2002a) and Ross, Fabian & Ballan- 
tyne (2002) suggested that they may strengthen the emission 
features and thus produce an apparent large reflection fraction. 
But several unknowns remain: since the reflection features are 
formed within 10 Thomson depths (tt) of the surface, what 
is the magnitude of any density change and where must it lie 
in the illuminated gas in order to alter the observed reflection 
spectrum? Do the accretion disk simulations produce the nec- 
essary fluctuations? The answers to these questions are crucial 
to the interpretation of X-ray spectra of accreting black holes. 

In this work, we investigate the properties of X-ray reflec- 
tion spectra from heterogeneous atmospheres. This paper (Pa- 
per 1) first considers a toy model where a density cut of vari- 
able depth and width has been sliced out of a constant density 
slab (§ 2). This will allow a systematic investigation of the ef- 
fects of a simple density inhomogeneity. We will then consider 
a more complicated structure produced by calculations of the 
photon-bubble instability in a radiation-dominated atmosphere 
with no turbulence (§ 3). Reflection spectra calculated from all 
these structures will then be compared with the standard con- 
stant density ones. In § 4 we discuss the results from this paper 
and form our conclusions in § 5. In a companion paper (Paper 
2), we apply the knowledge gained here to reflection by simu- 
lated radiation dominated accretion disks. 

2. TOY MODEL OF INHOMOGENEOUS REFLECTION 

The simplest test to see if density inhomogeneities affect re- 
flection spectra is to insert a discontinuous jump or drop in den- 
sity into an otherwise constant density model. In this section, 
we experiment with placing such steps with different widths at 
various places in a uniform slab. 

2.1. Model Setup and Assumptions 

The calculations used the reflection code described by Ross 
& Fabian (1993) and Ross et al. (1999), and the reader is re- 
ferred to those papers for the details on the computational pro- 
cedures. The models are computed by solving the coupled 
equations of radiative transfer and thermal and ionisation bal- 
ance in one spatial dimension, using the two-stream approxima- 
tion for the incoming radiation and a Fokker-Planck/diffusion 
method for the emergent spectrum. The solar abundances of 
Morrison & McCammon (1983) were employed in all the mod- 
els, as was an incident power-law radiation field (defined from 
1 eV to 100 keV) with photon index V = 2. The following ions 
were treated: C V - VII, N VI - VIII, O V - IX, Mg IX - XIII, 
S XI - XV and Fe XVI - XXVII. 

Each model started with a constant density slab (rin = 
10 15 cm~ 3 ) defined to have a total Thomson depth of tt = 6.1. 
At this optical depth almost all of the high energy photons will 
interact at least once with the gas. Beginning at a depth T ste p 
the density was then multiplied by a factor f p to create an ar- 
tificial drop or jump in the density profile at that position (tt 
remained fixed at 6.1 by adjusting the physical width of the 
altered zones). The increase or decrease in density was done 
discontinuously: there are no zones with densities intermediate 
between and f p nn. The width of the step was arbitrarily 
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defined as At = fcr stcp , where k is a constant. This results in 
the steps becoming wider as they are placed deeper which po- 
tentially allows for stronger observable effects to the reflection 
spectra. 

We chose to calculate reflection spectra for 4 different ioniza- 
tion parameters: £ =125, 250, 500, & 1000 erg cm s _1 . These 
values bracket the transition from a neutral 6.4 keV Fe Ka line 
to an ionized He-like 6.7 keV line (Ross et al. 1999). The soft 
X-ray features in the reflected emission, such as O VIII Lya, 
also vary greatly over this range of £ (Ballantyne, Ross & 
Fabian 2002). For each value of £, reflection spectra were cal- 
culated for steps placed at r stop =0.1, 0.5, 1, 1.5 & 2, with f p 
ranging from 10~ 3 , 10~ 2 , 0.1, 0.5, 2 & 10. The width of each 
step was controlled by the value of k, and for each (r s t ep , f p ) 
pair, models were computed for k = 0.2, 0.5, 1 & 2. Thus, 120 
different types of density steps were calculated for each £. All 
but two models converged successfully, and the resulting re- 
flection spectra were compared to an unaltered constant density 
model to assess the significance of any changes (§2.2). 

As an illustration of the technique, Figure 1 plots the tem- 
perature and density profiles as a function of tt for eight 
£ = 250 erg cm s _1 models. The results are shown for f p = 0.1 
in Fig. 1(a) and f p = 10.0 in Fig. 1(b). In each plot the dotted 
lines show the density profiles for four different steps (occur- 
ring at T stcp = 0.1, 0.5, 1 & 2) inserted into the slab, while the 
solid lines denote the final temperature profile. All the steps 
shown here have a width of At = 0.5T step (i.e., k — 0.5). The 
gas temperature reacts as expected: it increases by a factor of 
few within a density drop, and decreases within a jump. How- 
ever, the effect of an increase in nn seems to have a greater 
impact on the overall temperature profile than a sudden drop in 
density. Fig. 1(a) shows that outside of the step, the tempera- 
ture profiles are roughly similar, and indeed closely resembles 
the one from a constant density model (not shown). On the 
other hand, Fig. 1(b) shows that the temperature profiles are 
quite different for the four models. The differences between 
these two cases illustrates the major difference between an in- 
crease or a decrease in density. If the density changes by a fac- 
tor f p then the ionization parameter in that region is changed 
by 1/ f p . When f p < 1, the gas in the altered zones is more 
easily ionized which raises the temperature, but, most impor- 
tantly, the recombination and free-free emissivity rates (which 
scale as r? e w n^) are reduced. It is this increase in cooling 
when the gas becomes denser which has the largest impact on 
the temperature profiles seen in Fig. 1(b). However, the ionized 
gas in a low density region may be a strong emitter in some 
energy bands (e.g., the 6.7 keV Fe Ka line). Furthermore, a de- 
crease in density will also lower the opacity to radiation emitted 
from below. Therefore, we now turn to see how the impact of a 
sudden density change impacts the reflection spectrum. 

2.2. Results 

As a first step in determining the effects of a step in density 
on the reflection spectrum, we simply overplot the predicted 
emission with that from a constant density slab under the same 
illumination conditions. Four examples of this comparison, one 
for each value of £, are shown in the upper panels of Figure 2. 
Since in some cases any differences between the two spectra 
are difficult to see by eye, we plot in the lower panels the ra- 
tio between the constant density and the inhomogeneous model 
spectra. Prior to the ratio being calculated the T — 2 power- 
law was added to each spectrum so that each had a reflection 



fraction of unity. This was done in order for the plot to more 
accurately represent the ratio of two observed spectra (in AGN 

6 GBHCs the power-law continuum is observed along with the 
reflected emission). 

Figure 2, while only showing four examples, does demon- 
strate a number of useful points. First, although there are 
changes, the reflection spectra from the toy inhomogeneous 
models are not dramatically different from the constant density 
results; that is, they still look qualitatively like the spectra first 
presented by Ross & Fabian (1993). Secondly, the three vari- 
ables regarding the density step (its placement, depth/height, 
and width), as well as the value of £, all affect the deviations 
observed in the spectra. For example, the step illustrated in 
the £ = 250 erg cm s _1 panel of Fig. 2, which corresponds 
to one of the models shown in Fig. 1(a), results in a negligi- 
ble change to the reflection spectra above ^0.7 keV. The only 
major differences in the reflection spectrum are a decrease in 
intensity of low-energy Si XI, O V, and C V lines. In this case, 
the drop in energy was placed too deep into the layer for it to 
make much difference to the reprocessed emission. However, 
at larger ionization parameters, the incoming X-rays penetrate 
further into the slab, and so a step at T stop = 1 can have a 
large impact on the reflection spectrum. This is seen in the 
£ = 500 erg cm s" 1 panel of Fig. 2, but more dramatically 
in the £ = 1000 erg cm s" 1 plot. In both these cases, the 
low density step alters the reflection spectra so that they have 
a shape that resembles ones with a larger value of £ (less ab- 
sorption, broader spectral features). These two plots show that 
the reflection spectra from a inhomogeneous slab can be a mix- 
ture of ionization parameters, similar to ones calculated from 
hydrostatic atmospheres (Ballantyne et al. 2001b). Of course, 
when the step causes an increase in density, as is shown in the 
£ = 125 erg cm s _1 panel, then a spectrum typical of a lower 
ionization parameter (stronger absorption and low-energy emis- 
sion lines) is mixed in with the results. 

Although Fig. 2 does show that density inhomogeneities can 
cause observable effects in the reflection spectra, the results are 
largely qualitative and selective (there are ~ 475 other possible 
panels). To better quantify the results, and to obtain a more 
general overview of the parameter space, we calculated the 
Fe Ka and O VIII Lya EWs from the total (reflected+incident) 
spectrum for each two-density model and compared it to the 
values from the corresponding constant density models. The 
EWs were calculated by integrating the spectrum between 6 and 

7 keV for Fe Ka, and between 0.6 and 0.7 keV for O VIII Lya. 
Figure 3 shows the results of this exercise, where for each £ we 
plot the EW ratio (defined as two-density model/constant den- 
sity) for both the Fe Ka (black) and O VIII (blue) lines as a 
function of T stcp . Each plot contains 6 panels for the different 
values of f p , and the different style of lines within the panel de- 
notes a particular value of k. There is plenty of information in 
this figure, and it is worth going through it in detail, but before 
doing so there is one general point that can be made. In almost 
every instance the magnitude of the EW ratio shown in Fig. 3 
is proportional to the value of k; that is, the wider the step, the 
greater the change to the reflection spectrum for those values of 
f p and T stop . 

Starting with the £ = 125 erg cm s _1 plot, we find that be- 
cause of the relatively weak illumination, the density change 
generally only affects the EWs for T stop < 0.5. The largest 
changes to the Fe Ka EW are when f p = 0.1, 2 or 10. In 
the first case, the increase in EW is due to the appearance of 
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a 6.7 keV component to the Fe Ka line, while the later two 
instances are a result of a stronger 6.4 keV component. The 
reason why the ionized line appears in the f p = 0.1 models 
but not in the f p = 0.01 or 0.001 cases is that the Fe is al- 
most entirely fully stripped in the very low density zones and 
so does not contribute much emission (as do the other metals). 
The O VIII EW varies only slightly, except for when f p = 10 
where it drops by <~ 30% due to a low ionization fraction in the 
overdense region. 

A similar pattern is seen in the £ = 250 erg cm s _1 models, 
with the changes now occurring for r stcp < 1. There are sig- 
nificant variations in the Fe Ka EW for most values of f p , with 
the largest changes (> factor of 8) occurring in the f p = 0.1 
case. The reflection spectrum for the £ — 250 erg cm s _1 con- 
stant density model (seen as the dashed line in Fig. 2) has a 
very weak Fe Ka line due to Auger destruction (Ross, Fabian 
& Brandt 1996). Therefore, changes to the density which re- 
sult in either a 6.7 keV line (f p = 0.1) or a strong 6.4 keV line 
(f p = 10) cause a significant increase in the EW. At this value 
of £ it seems to be difficult to increase the O VIII EW by den- 
sity steps, but it can be lowered by nearly a factor of two when 

f P = io. 

Increasing the illumination level so that £ = 500 erg cm s 1 
results in some changes in the reflected emission with a step 
down at T stcp ~ 2, but, from Fig. 3, we see that the most sub- 
stantial changes still occur at T stcp < 0.5. The Fe Ka EW 
differs significantly from the constant density value only when 
f p = 0.1 or 0.5. In these cases the iron in the low density zones 
becomes highly ionized, but not fully stripped, increasing the 
column of the He-like species. The line is also broader due 
to Compton down-scattering. Again, we see that the spectrum 
takes on the properties of one from a larger £. This is also re- 
flected in the changes in the O VIII EW, which is much smaller 
than the constant density model for f p = 0.001-0.5. How- 
ever, the opposite behavior (a larger O VIII EW and a smaller 
Fe Ka EW) is seen when the step causes an increase in den- 
sity and low-£ features are mixed into the reflection spectrum. 
In most of those models, the 6.7 keV Fe line and the O VIII 
Ly a line are reduced and enhanced in intensity, respectively, 
although an interesting counter-example is seen when f p = 10 
and r stop = 0.5. In this model, the overdense region is placed 
at a depth so that both a 6.4 keV and a 6.7 keV Fe Ka line is 
emitted, causing an increase in the EW. Similarly, the O VIII 
EW falls because of a significant transition to O VII emission. 

Finally, the most highly ionized runs (£ = 1000 erg cm s _1 ) 
show little variation in the Fe Ka EW when f p < 1. The 
only exceptions are when the density deficit is near the surface 
and the high ionization parameter in these zones decreases the 
amount of 6.7 keV emitting ions. As a result, the line EW falls 
by about 25%. On the other hand, when f p = 10, the Fe Ka 
EW can drop by over a factor of two because the overdense re- 
gion stops the ionization front at a much shallower depth in the 
slab. The O VIII EW is always smaller than its constant density 
value when f p < 1, but can increase by over a factor of 2 (up to 
110 eV) when there is a density enhancement just beneath the 
surface. 

In summary, we find that toy models which place a top-hat 
function of variable width and depth into a constant density 
slab, can produce observable effects on the resulting reflection 
spectra. Generally the largest impact on the spectrum, as mea- 
sured by the Fe Ka and O VIII EWs, occurs when the step is 
placed at T stcp ~ 0.5 into the layer. Furthermore, changing 



the density by only relatively small factors of 2-10 can produce 
significant effects. The resulting reflection spectra are reminis- 
cent of the one produced by hydrostatic models in so far as they 
show features from a mixture of ionization parameters. 

2.3. Fits with Constant Density Reflection Spectra 

The EW ratios shown in Fig. 3 indicate that density inhomo- 
geneities in the surfaces of accretion disks may have an observ- 
able effect on their reflection spectra. Those measurements, 
however, only use a fraction of the information contained in 
the model spectra. In this section, we take several example 
reflection spectra calculated from the two-density models and 
fit them with ones from the constant density models. This 
method will allow us to quantify any differences between the 
two scenarios by using all the information over a given energy 
range. The exercise may also be indicative of how real inhomo- 
geneities will manifest themselves when fitting data with simple 
models. 

We chose models that showed large changes in its Fe Ka or 
O VIII EW from each of the four ionization parameters con- 
sidered (the exact model parameters are listed in the note to 
Table 1). Simulated observations were made with these model 
spectra by using the 'fakeit' command in XSPECv.l 1.2.0bp (Ar- 
naud 1996). The canned XMM-Newton EPIC-pn response ma- 
trix epn_sw2 0_sdY9 . rmf was used to generate the data as- 
suming an exposure time of 40 ks. As with the EW calculation 
in the previous section, the incident power-law was added to the 
reflection spectra so that each simulated spectrum had a reflec- 
tion fraction, R, of unity. Each faked dataset was then fit with a 
grid of constant density models (calculated with the same ver- 
sion of the reflection code) between 0.2 and 12 keV as well as 
between 2 and 10 keV. The fit parameters were the photon in- 
dex T, the ionization parameter £, the reflection fraction R, and 
an arbitrary normalization. The results are listed in Table 1, and 
the residuals are shown in Figure 4. 

This particular £ = 125 erg cm s _1 model was chosen be- 
cause the Fe Ka EW increased by over a factor of 2 due to the 
appearance of a 6.7 keV line. In the 0.2-12 keV energy band, 
the constant density model does a decent job of recovering the 
ionization parameter, but the value of R is too low. In con- 
trast, the best fit in the 2-10 keV range has a much larger £ 
(466 erg cm s _1 ) and R. Evidently, the spectrum in this energy 
band mimics one from a more ionized reflector. The residuals 
to the 0.2-12 keV fit are shown in Fig. 4 and shows that the 
low value of R is likely due to the overprediction of the O VIII 
line (as seen in Fig. 3). Moreover, the ionized Fe Ka line can- 
not be accounted for by this low £. Ignoring the low energy 
data increases the reflection fraction as the spectrum attempts 
to account for the enhanced Fe emission. In both cases, the 
continuum is fit adequately, but the differences in line emission 
produce significant residuals. 

Statistically, the fit to the £ = 250 erg cms -1 two-density 
model over the 0.2-12 keV band is much worse than for the pre- 
vious case. The best fit ionization parameter is lower than the 
true one, and so overpredicts the strength of the soft emission 
lines between 0.3 and 0.6 keV (Fig. 4). According to Fig. 3, the 
Fe Ka EW in this two-density model was ~2 times greater than 
the corresponding constant density model (although at this £ the 
EW is only 20 eV), but that the change in the O VIII EW was 
marginal. Thus, in this case we find that the overdense region 
lowers the effective ionization parameter of the reflection spec- 
trum, but this is only discernible through a fit to the soft X-ray 
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emission. Indeed, we find a very good and statistically accept- 
able fit to this two-density model over the 2-10 keV band. The 
best fit £ and R are only slightly larger than their true values to 
account for the enhanced Fe Ka line. 

The two-density model chosen from the £ = 500 erg cm s _1 
suite of calculations is another overdense model with the step 
appearing at T stcp = 0.5 below the surface. As discussed in 
§ 2.2, this model includes both 6.4 and 6.7 keV iron emission 
lines, as well as a weaker-than-expected O VIII line. Interest- 
ingly, this model gives the worst result when fit by the con- 
stant density models, with the reduced \ 2 > 3. As in the 
£ = 250 erg cm s _1 model, the poor fit is caused almost entirely 
by an overprediction of the soft X-ray spectral features. The 
majority of the 0.2-12 keV continuum is well fit by an ionized 
reflector (the best fit £ is ^630 erg cm s _1 ), but the emission 
lines cannot be accounted for by this model, resulting in the low 
value of R. The overdense region sitting just below the surface 
of the reflector produces features that cannot be fit by a single-^ 
model. The constant density models do a much better job in the 
2-10 keV band, where the best fit ionization parameter is nearly 
exactly the true one. However, the reflection fraction must be 
lowered to fit the weaker 6.7 keV Fe line. There are still some 
residuals in this fit between 7 and 8 keV where the overdense 
region in the two-density model causes additional absorption 
than what is predicted by the constant density model. 

Our final example considers the case where higher-^ fea- 
tures are mixed into the reflection spectrum. The £ = 
1000 erg cm s _1 model used for fitting had a step with a density 
drop of f p = 0.001 at r stop = 0.5. The lower density mate- 
rial became highly ionized and weakened both the Fe Ka and 
O VIII lines (Fig. 3). The two constant density fits both show 
the effects of this dilution, as the best fit £ is greater than the 
true value and R < 1. The residuals shown in Fig. 4 show that 
an overprediction of soft X-ray features is again the cause for 
the poor fit. There is also a significant difference in the curva- 
ture of the two spectra noticeable between 1 and 3 keV. This is 
caused by the lower absorption in the two-density model which 
flattens the continuum. 

3. PHOTON BUBBLES AND REFLECTION 

The idealised two-density models discussed in the previ- 
ous section show that inhomogeneities in accretion disk atmo- 
spheres may affect the reflection spectrum. We next exam- 
ine reflection from density structures arising in hydrodynamical 
disk models. Effects of the ionization processes on the density 
structure are neglected. Nevertheless, the density arrangement 
is completely novel in terms of reflection calculations, as it is 
the first to be considered which did not arise from an ad hoc a 
priori assumption (such as constant density or hydrostatic bal- 
ance). 

The two-density results indicate that the reflection spectrum 
is affected by inhomogeneities within the outermost few Thom- 
son depths of the disk atmosphere. By contrast, fluctuations re- 
sulting from the MRI may be strongest on scales comparable to 
the disk thickness (Hawley, Gammie & Balbus 1995; Stone et 
al. 1996). The effects of photon bubbles are therefore consid- 
ered instead. The PBI grows most rapidly at short wavelengths 
(Blaes & Socrates 2001), and in Shakura-Sunyaev models typ- 
ically grows fastest in the disk surface layers (N.J. Turner et al. 
2004, in preparation). Exploring X-ray reprocessing from pho- 
ton bubbles is likely to be a good first step toward understanding 
the reflection properties of radiation-dominated disks. 



The photon bubble instability was first discussed by Arons 
(1992) in the context of accreting magnetized neutron stars, and 
was shown by Gammie (1998) to occur in the inner, radiation- 
supported regions of magnetized accretion disks. For magnetic 
pressures exceeding the gas pressure, and perturbation wave- 
lengths shorter than the gas scale height, the instability oper- 
ates when disturbances in radiative flux displace gas along field 
lines, leading to propagating density variations which grow 
over time. Linear growth is fastest at wavelengths shorter than 
the gas pressure scale height, but is absent at wavelengths with 
optical depths much less than unity, as these support no flux 
variations. An approximate criterion for instability is that pho- 
ton diffusion carry a greater energy flux than radiation advected 
at the gas sound speed (Blaes & Socrates 2003). 

Reflection spectra were calculated using a density profile 
taken from a two-dimensional radiation-MHD simulation of 
the growth of photon bubbles in a Shakura-Sunyaev model 
disk. The simulation was carried out with the Zeus MHD code 
(Stone & Norman 1992a,b) and its flux-limited radiation dif- 
fusion module (Turner & Stone 2001). Opacities due to elec- 
tron scattering and bremsstrahlung were included. The domain 
was a small patch of the disk, centered 34 i?s from a black 
hole of 10 8 M Q . The patch extended 1.08 Shakura-Sunyaev 
scale-heights either side of the midplane, and its width was one- 
quarter of its height. The grid consisted of 64 x 256 zones. The 
initial density and temperature profiles were those of a Shakura- 
Sunyaev model with a — 0.06, and accretion rate 12% of the 
Eddington rate for a 10% luminous efficiency. During the sim- 
ulation, no a-viscosity was applied. Differential rotation was 
neglected, so that the MRI was absent, and the material was 
allowed to cool by radiative losses. Photons initially diffused 
from midplane to boundary in nine orbits. The calculation was 
started with a magnetic field having a pressure 10% of the mid- 
plane radiation pressure, and inclined at 45°. After 3.4 orbits, 
the PBI had developed into trains of propagating shocks, with 
density contrasts of up to two orders of magnitude. The verti- 
cal density profile used in the subsequent reflection calculations 
was taken from an arbitrary vertical ray, and is shown in the in- 
sets to Figure 5. 

3.1. Results 

As with the two-density models, this new density struc- 
ture was then irradiated by a power-law (T = 2) of X-rays. 
In this case, there is no unique ionization parameter so in- 
stead the models are distinguished by the incident X-ray flux: 
F x = 10 13 , 2.5 x 10 13 , 5 x 10 13 , and 10 14 erg cm~ 2 s" 1 . 
The Shakura-Sunyaev flux for the assumed disk parameters is 
3 x 10 13 erg cm~ 2 s _1 . All four models converged, and the 
resulting reflection spectra are presented in Figure 5. The gas 
temperature and number density are plotted in the insets to each 
panel, and illustrates how the temperature structure for each 
case was affected by the density profile. 

The latter three spectra all exhibit strong 6.7 keV emission 
lines from He-like iron, while, unusually, the most weakly illu- 
minated model produces both a 6.4 and a 6.7 keV line of about 
equal strength. The rapid drop in density below one Thomson 
depth causes the surface layers to be easily ionized and able to 
produce a 6.7 keV emission line even with a relatively small in- 
cident flux. On the other hand, the overdense region at tt ~ 1 
produces a lot of soft X-ray line emission, even for highly illu- 
minated situations. Thus, the emission below 2 keV is consis- 
tently important, despite the incident flux varying by an order 
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of magnitude. 

3.2. Fits with Constant Density Reflection Spectra 

To quantify any effects the photon bubbles may have on the 
reflection spectra we fit the models with a grid of constant 
density spectra. The procedure was exactly the same as with 
the two-density models in § 2.3: the power-law was added 
to the photon bubble spectra and then the sum was used as 
the basis for a 40 ks XMM-Newton simulation. As before, 
constant density fits were performed in both the 0.2-12 keV 
and 2-10 keV energy bands. A constant density grid with 
riH = 1.4 x 10 12 cm~ 3 was used for fitting the photon bubble 
reflection spectra in order to offset any changes in the reflection 
spectra due solely to the lower density. The best fit parameters 
for each value of Fx are listed in Table 2, and the residuals are 
shown in Figure 6. 

The constant density models had difficulty fitting the spectra 
over the wide 0.2-12 keV energy range, with most of the trou- 
ble arising from the soft X-ray lines. The best \ 2 was found 
with the Fx = 10 13 erg cm~ 2 s _1 model, where one ionization 
parameter seemed to provide a good fit to the spectrum at ener- 
gies > 0.7 keV. The best fit £ was indicative of an ionized slab, 
and, indeed, the model could not account for the small 6.4 keV 
line predicted by the spectrum. Despite this problem, the fit was 
able to recover the correct values of both V and R. The fit had 
difficulty below 0.5 keV, however, where it overpredicted the 
emission. This is most likely a result of extra absorption (over 
that predicted for the fitted £) due to the density enhancement 
just beneath the surface in the photon bubble model. 

Fits to the remaining three models in this energy band all re- 
sulted in reduced x 2 > 2. A single ionization parameter could 
not simultaneously account for the shape of the continuum, and 
the strong soft X-ray and 6.7 keV Fe Ka lines. The best fit val- 
ues of R were always less than the 'true' value of unity. In these 
cases the photon bubble models predicted significant emission 
from a variety of ionization parameters, so that the resulting 
spectrum exhibited a mixture of features. As a result, a single 
ionization parameter fit was not adequate. 

The fits to the photon bubble models improved markedly in 
the 2-10 keV energy band with £ increasing substantially in 
order to fit the Fe Ka line. As a result, all four fits were sta- 
tistically acceptable. Except for the Fx = 10 13 erg cm~ 2 s _1 
model, the reflection fractions again are still underestimated by 
the constant density fits with the best fit value of R decreasing 
with the Fx of the model. A greater than expected amount of 
absorption is the likely explanation for this effect. The den- 
sity enhancement at tt ~ 1 in the photon bubble model causes 
more continuum absorption from oxygen then what the reflec- 
tion models expect (given, e.g., the highly ionized Fe Ka line). 
This enhancement becomes more important as Fx is increased 
and spectral features are formed both above and below it. The 
extra absorption decreases the magnitude of the continuum and 
emission, and thus results in a lower R when fit with the con- 
stant density models. Aside from the low value of R, reflection 
from a uniform slab does a good job describing the photon bub- 
ble spectra between 2 and 10 keV. 

4. DISCUSSION 

Currently, the interpretation of the reflection continuum ob- 
served in the X-ray spectra of AGN is the primary means to 
glean information regarding the accretion geometry and its ra- 
diative environment. It is therefore vital that the models of re- 



flection, on which such conclusions are based, explore a num- 
ber of different physical situations which may be relevant to 
accretion physics. Motivated by the inhomogeneous nature of 
recent numerical simulations of radiation-dominated accretion 
disks, we have begun a study of examining the consequences of 
such density changes to the X-ray reflection spectrum. 

Our first results, presented in the previous sections, consid- 
ered reprocessing from a two-density model, where density 
steps or jumps were inserted into a uniform slab, and a slice 
from a PBI simulation. The new reflection spectra were then 
compared with ones calculated assuming a constant density at- 
mosphere, as these are most frequently used in data analysis. 
We found that the density inhomogeneities can result in ob- 
servationally important differences between the two cases. As 
in models with a hydrostatic atmosphere, the reflection spectra 
can no longer be described by a single ionization parameter, but 
exhibit features from a mixture of ionization states. However, 
while in the hydrostatic case the spectra can be adequately de- 
scribed as a diluted constant density model because of a diffuse 
hot scattering layer on the surface (Ballantyne et al. 2001b), the 
inhomogeneous models can show emission from both low and 
high £ material simultaneously, depending on the nature of the 
density change. This gives the spectra a complexity and rich- 
ness that the hydrostatic models lacked. 

The magnitude and characterization of the change in the re- 
flection spectrum depends on both the nature of the inhomo- 
geneity and the strength of the incident X-rays. For the two- 
density models, where the calculations could be distinguished 
by the original ionization parameter £ of the uniform slab, the 
greatest impact on the Fe Ka and O VIII Lya EWs occurred 
whenT s tcp < 1 for£ < 1000. For larger ionization parameters, 
steps deeper into the gas could make an impact. This correla- 
tion is not expected to continue indefinitely however. For ex- 
ample, if an atmosphere had a drop in density at tt ~ 3 and it 
was illuminated to the extent that gas was ionized down to this 
depth, then the additional line emission or absorption caused by 
the inhomogeneity would be smeared out by Compton scatter- 
ing while escaping the layer. Thus, a general conclusion seems 
to be that density inhomogeneities must be within 2 Thomson 
depths of the surface to have any impact on the reflection spec- 
trum. 

However, there must also be a lower-limit to the depth of 
any density change in the reflecting medium for it to be effec- 
tive in altering the spectrum. In the extreme case of cold re- 
flection (e.g., £ w 10) ionization effects are important only at 
tt < 0.1 from the surface. A density inhomogeniety at such a 
depth would change such a small amount of gas that it would 
have a negligible effect on the resulting spectrum. 

The results of the two-density models also indicate that rela- 
tively small changes in density (say, Sp/p ~few) can alter the 
reprocessed emission. Again, this depends on the effective ion- 
ization parameter of the layer, but, as was seen in § 3, under- 
dense regions beneath the surface of a moderately irradiated 
atmosphere can result in a substantial 6.7 keV iron line. But, 
if that region was underdense but a large amount, such as 100- 
1000 times lower, then it will be ionized to such an extent that 
it has very little effect on the outgoing spectrum. Similarly, an 
overdense region below the surface may enhance the 6.4 keV 
Fe Ka line. Analogous arguments also apply to the soft X-ray 
emission lines such as O VIII. 

How realistic are these requirements on the density inhomo- 
geneities? It certainly appears plausible, even likely, that accre- 
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tion flows will naturally generate density contrasts greater than 
a few, especially in the radiation dominated regime where the 
photon bubble instability can enhance already existing dumpi- 
ness. What remains unknown is if the heterogeneous nature 
of the flow exists to small enough scales that sufficient inho- 
mogeneity remains within a couple of tt from the disk photo- 
sphere. Only very high resolution simulations of this region of 
accretion disks will be able to definitively answer that question. 

Assuming for the moment that the inhomogeneities dis- 
cussed in this paper occur in reality, what are the immediate 
observational consequences? Perhaps the most interesting re- 
sult is that the density jumps can cause the Fe Ka line to vary 
in a way that is completely disconnected from the X-ray con- 
tinuum. One of the most well known puzzling properties of 
the possible Fe Ka 'disklines' is their lack of response to the 
variable continuum (e.g. Chiang et al. 2000). The line flux 
typically does vary, but it is often not directly correlated with 
the continuum (Iwasawa et al. 1996; Weaver et al. 2001; Wang 
et al. 2001; Vaughan & Edelson 2001; Markowitz et al. 2003; 
Iwasawa et al. 2003). One very striking example, which was 
presented by Petrucci et al. (2002), is an XMM-Newton ob- 
servation of Mkn 841 that was split into two parts, separated 
by about 15 hours. These authors found that the EW of the 
narrow Fe Ka dropped by ~ 2 between the two observations 
while the continuum changed by ~ 10-20%. Although, we are 
not explicitly modeling this source, Fig. 3 shows that density 
steps in the surface can give rise to exactly such changes in the 
Fe Ka EW, completely independently of the X-ray continuum. 
A clear way to test if this is the correct model would be to ex- 
amine the changes of other emission features in the spectrum, 
but this would require much more sensitive data. The combina- 
tion of density inhomogeneities with the non-monotonic evo- 
lution of the Fe Ka EW due to ionization effects (Ballantyne 
& Ross 2002) may naturally result in a poor correlation be- 
tween the observed flux and the Fe Ka EW. However, since the 
observed X-rays must be averaged for many kiloseconds be- 
fore spectral analysis, which is longer than the timescale for a 
change in the density inhomogeneities, as well as the the ion- 
ization/recombination timescale, a relationship between the line 
and continuum may still be uncovered, but at a weaker level 
than what may have been previously expected. 

If the density inhomogeneities beneath the photosphere are 
relatively small in scale, which should be true for thin disks 
(where H < R, and H is the disk scale height and maxi- 
mum size of any density fluctuations), then any resulting en- 
hancement or change to the Fe Ka line would only occur to 
a narrow part of the overall line profile, which is determined 
by the illuminating emissivity. Thus, the appearance of narrow 
Fe Ka lines, such as those recently inferred from observations 
of NGC 3516 (Turner et al. 2002), NGC 7314 (Yaqoob et al. 
2003), and Mrk 766 (Turner, Kraemer & Reeves 2003) could 
conceivably originate from density inhomogeneities in the disk. 

Variable Fe Ka lines from turbulent accretion disks have also 
been considered by Armitage & Reynolds (2003). These au- 
thors assumed the line emissivity was proportional to the local 
integrated magnetic stress in their numerical simulations, as op- 
posed to our method of assuming a constant illumination and a 
variable density structure. In reality, a mixture of the two effects 
would be expected to be ongoing, which should be investigated 
in future work. 

The strength of the soft X-ray emission lines are also affected 
by density inhomogeneities. These lines are one of the domi- 



nant sources of cooling for the X-ray heated gas when it reached 
< 10 6 K. Thus, if overdense or underdense regions are intro- 
duced into the layer it can alter the rapidity at which the gas 
cools, thereby changing the emission features. Interestingly, 
these changes may occur only in the soft band, with little ef- 
fect at higher energies. Soft X-ray lines, especially the Lya 
lines of O VIII, N VII and C VI, have received some attention 
recently with the claims that they have been observed to be rela- 
tivistically broadened in the gratings spectra of some Seyfert Is 
(Branduardi-Raymont et al. 2001; Mason et al. 2003; Sako et 
al. 2003). The exact strength of the lines have been the subject 
to some debate in the literature (Lee et al. 2001; Ballantyne et 
al. 2002; A.K. Turner et al. 2003), with the difficulties arising 
because complex warm absorption features must be taken into 
account in the spectral modeling. The results from the two- 
density models show that density inhomogeneities can both en- 
hance and diminish the EWs of the soft X-ray lines. 

The constant density fits to both the two-density models and 
the spectra computed from the photon bubble structure resulted 
mostly in reflection fractions much less than the true value of 
unity. If one interprets R as a measure of the solid angle sub- 
tended by the disk as seen from the X-ray source, then R < 1 
would indicate some form of truncated accretion flow (e.g., 
Zycki, Done & Smith 1997, 1998, 1999; Done & Zycki 1999; 
Eracleous, Sambruna & Mushotzky 2000; Gliozzi, Sambruna 
& Eracleous 2003). Correctly interpreting a fitted value of R 
is fraught with difficulty since fitting neutral reflection models 
to ionized accretion disks will also give R < 1 even if R = 1 
(Ballantyne et al. 2001b; Done & Nayakshin 2001). Here, we 
have shown that ionized disk models can produce erroneous R 
values when there are density inhomogeneities. It is therefore 
important that other arguments be used (e.g., Barrio, Done & 
Nayakshin 2003) before drawing conclusions on the accretion 
geometry from a fitted value of R. 

5. CONCLUSIONS 

Our conclusions from this first paper in our study on reflec- 
tion from heterogeneous accretion disks can be summarized as 
follows: 

1 . The reflection spectrum from accretion disks can be al- 
tered by clumps or voids in the gas just beneath the pho- 
tosphere. The effects range from negligible to large and 
depend on the structure of the inhomogeneity and the 
illuminating continuum. 

2. The greatest effects on the reflected emission for most 
ionization parameters occurs when the density change 
is within ~ 2 Thomson depth of the surface. 

3. The change in density does not have to be large, even an 
increase or a decrease by a factor of a few can signifi- 
cantly alter the spectrum from that of a constant density 
model. 

4. Density inhomogeneities beneath the disk surface are a 
possible explanation for the apparent disconnectedness 
between the Fe Ka line and the X-ray continuum. 

5. The soft X-ray emission lines (below 1 keV) are more 
sensitive to the presence of an inhomogeneity than the 
harder emission. 

6. Fitting the spectra produced from models with inho- 
mogeneities with constant density models shows that a 
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mixture of ionization parameters are present, which of- 
ten results in reflection fractions smaller than unity. 



7. Constant density models may still be a good means of 
parameterizing reflection spectra, but only for energies 
> 2 keV. They will have much more difficulty with 
broadband spectra that cover the soft X-ray lines emit- 
ted from the disk. However, in practice this may only be 
a problem for sources with weak or non-existent warm 
absorbers. 



In Paper 2, we continue our investigation by computing re- 
flection spectra from the predictions of 3-D radiation dominated 
accretion disk simulations. Future work will also include a con- 
sideration of the dynamical effects of the incident X-rays and 
multi-dimensional reflection. 
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Table 1 



Results of fitting a sample of two-density models with constant density reflection spectra. 



£modcl 


N 


r 


0.2- 


12 keV 
R 


X 2 /d.o.f. 


r 


2.0-10 keV 
C R 


x 2 /d.o.f. 


125 
250 
500 
1000 


jQ-23 
1Q -23 

5 x 10- 24 
2 x 10~ 24 


9 fy-,, +0.004 
z " uz ' i -0.003 

1 mfi+ 002 
z.uzo_ 004 

2.0 

2 000+ 004 


147±3 
190±3 

631 
1244+1? 


o.79tg;S 

0.92t°;g 2 

0.654 
0.75±0.02 


2672/2358 
4265/2358 
7833/2358 
2994/2358 


2.03±0.01 
2.02±0.01 
2.03±0.01 


46611 
327+1° 
507+ 2 ° 
128211™ 


1.2±0.2 
1.1±0.1 

0.88l°;° 8 9 
0.81±0.05 


1684/1595 
1590/1595 
1638/1595 
1639/1595 



Note. — The step parameters for the two-density model spectra are (a) £ = 125, f p = 0.1, r stC p = 0.1, k = 2; (b) £ = 250, f p = 10, r s to P = 0.5, 
k = 0.2; (c) £ = 500, f p = 10, r st0 p = 0.5, k = 1; (d) £ = 1000, / p = 0.001, r stC p = 0.5, fc = 2. JV is the normalization of the two-density models 
in XSPEC. R is the reflection fraction defined as total=incident+i?x reflected. £ is in units of erg cm s _1 . The error-bars are the 2a uncertainty for the 
parameter of interest. The simulated data were constructed using the XMM-Newton response matrix epn_sw2 0_sdY9 . rmf and assumed an exposure 
time of 40 ks. 



Table 2 

Results of fitting the photon bubble reflection spectra with constant density models. 



Fx 


N 


r 


0.2- 


-12keV 
R 


X 2 /d.o.f. 


r 


2.0-10 keV 
C R 


x 2 /d.o.f. 


10 13 
2.5 x 10 13 
5 x 10 13 
10 14 


2 x 10~ 20 

1Q -20 

5 x ltr 21 

2 x 10~ 21 


1.9881°:°*° 

1.950 
1.970 
2.029 


56llg 

565 
593 
491 


1.01±0.01 
0.798 
0.926 
0.611 


4679/2358 
5821/2358 
6651/2358 
5995/2358 


2.04±0.01 
2.03±0.01 
2.00±0.01 

1 QR7+0.007 
i.yo l _Q.008 


60111 2 
1297tgf, 2 
1340lf 3 8 8 
2679t 334 


o.96l8:iJ 

0.73i»; 3 
0.52±0.03 
0.46±0.03 


1632/1595 
1695/1595 
1561/1595 
1627/1595 



Note. — iV is the normalization of the photon bubble models in XSPEC. R is the reflection fraction defined as total=incident+i?x reflected. £ is in units 
of erg cm s _1 . Fx is in units of erg cm~ s _1 . The error-bars are the 2a uncertainty for the parameter of interest. The simulated data were constructed 
using the XMM-Newton response matrix epn_sw2 0_sdY9 . rmf and assumed an exposure time of 40 ks. 
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FIG. 1 . — Gas temperature (solid lines) & number density (dotted lines) as a function of tt for eight § = 250 erg cm s _1 models, (a) Results for the case when 
f p =0.\. Steps were placed at T stC p = 0.1, 0.5, 1 & 2, each with At = 0.5r stC p. The gas temperature increases at the position of the drop in density, but otherwise 
there is little impact in the overall shape of the temperature profile, (b) As in (a), except f p = 10. In this case, the temperature falls at the jump in density, but, 
because of the increased cooling rates at lower T, it does not recover to a similar overall profile. 
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FIG. 2. — Four examples of how the reflection spectra computed from the toy models compare with those from unaltered constant density slabs. The top panel 
of each plot shows the reflection spectra as the solid line with the constant density one denoted by the long-dashed line. The short-dashed line shows the T = 2 
incident power-law used for both calculations. The lower panel plots the ratio of the two reflection spectra defined as (constant density/two-density model). With 
AGN the power-law is also observed along with the reflection spectrum, so the F = 2 spectrum was added to the reprocessed emission prior to the ratio being 
calculated (i.e., they both have a reflection fraction of unity). In this way, we are comparing the spectra as they may be observed. 
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FIG. 3. — These plots show as a function of r s top the ratio of the Fe K« (black) and O VIII (blue) Lya equivelant widths from the toy inhomogeneous models to 
those from the constant density models. The six panels in each plot descibe the results for each value of f p , while the linestyles denote the different values of k: 
solid (k = 0.2), dotted (k = 0.5), short-dashed (k = 1) and long-dashed (fe = 2). The EWs were calculated from the total (reflection+incident) spectra. See the 
text for the discussion. 
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FIG. 5. — X-ray reflection spectra (solid line) calculated from the outer 10 Thomson depths of a photon bubble simulation. The short-dashed line denotes the 
r = 2 power-law that was incident on the material. The insert shows the gas temperature (solid line) and number density (dashed line) for each model. Of 
course, the density is the same for each calculation. The different plots show the results for various incident X-ray fluxes: (a) Fx = 10 13 erg cm~ 2 s _1 , (b) 
F x = 2.5 X 10 13 ergcm" 2 s _1 . 
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FIG. 6. — Simulated count spectra and residuals (in units of standard deviations) of models computed using the density cut from the photon bubble simulation, 
and then fit with constant density reflection spectra. The fit parameters are given in Table 2. 



